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Abstract. In this paper we present results on transverse mass spectra and Hanbury- 
Brown and Twiss correlation functions of negatively charged hadrons, which are 
expected to be mostly 7r~, measured in Pb-Pb collisions at 40 A GeV/c beam 
momentum. Based on these data, the collision dynamics and the space-time extent of 
the system at the thermal freeze-out are studied over a centrality range corresponding 
to the most central 53% of the Pb-Pb inelastic cross section. Comparisons with 
freeze-out conditions of strange particles and HBT results from other experiments 
are discussed. 



PACS numbers: 12.38.Mh, 25.75.Nq, 25.75.Ld, 25.75.Dw 
1. Introduction 

Ultra-relativistic collisions between heavy ions are used to study the properties of 
nuclear matter at high energy density. In particular, lattice QCD calculations predict 
a transition from confined hadronic matter to a state of deconfined quarks and gluons 
at a critical energy density around 1 GeV/fm 3 pQ. For recent reviews of experimental 
results and theoretical developements see references [2]. 

The momentum distributions of the particles emerging from the Pb-Pb interactions 
are expected to be sensitive to the collision dynamics. In particular, collective dynamics 
in the transverse direction is of major interest since it can only arise by the buildup 
of a pressure gradient in that direction; this in turn would be strongly suggestive of 
thermal equilibration of the nuclear matter. Indeed such an effect has already been 
observed at the highest SPS and RHIC energies. The shapes of the m T = \/pj + m 2 
distributions are expected to be determined by an interplay between two effects: the 
thermal motion of the particles in the fireball and a pressure-driven radial flow, induced 
by the fireball expansion. In reference [3] we have analyzed the m T distributions of 
strange particles (A, E~, Q~, their anti-hyperons and Kg) based on the blast- wave 
model [1], a parameterized model inspired to hydro-dynamics. Due to the large number 
of particle species considered in that analysis, a simultaneous fit to the strange particle 
spectra allowed to disentangle the effect of the thermal motion from that of the collective 
expansion. With one particle species only, e.g. the negative pion, an oblongated 
confidence region can be obtained for the pair of freeze-out parameters temperature 
(T) and average transverse flow velocity ((P±)). 

A systematic study of the space-time extent and the dynamical behavior of the 
fireball at thermal freeze-out can be obtained via the identical particle interferometry 
technique, first used by Hanbury-Brown and Twiss (HBT) [5]. The width of 
the correlation peak at vanishing relative momenta reflects the so-called length of 
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homogeneity (also called the "HBT radius") of the particle emitting source. Only in 
static sources can the homogeneity length be interpreted as the true geometrical size of 
the system. In a dynamic system, the occurrence of space-momentum correlations of 
the emitted particles due to collective expansion generally leads to a reduction of the 
observed HBT radii. The degree of reduction depends on the gradients of the collective 
expansion velocity and on the thermal velocity of the particles at thermal freeze-out. 
A differential analysis of the HBT correlations in bins of the pair momentum (in the 
longitudinal and transverse directions) thus provides valuable information both on the 
spatial extent and on the properties of the collective expansion of the system. Combining 
single particle spectra and two particle HBT correlations allows us to disentangle the 
collective dynamics from the thermal motion relying on one particle species only. 

In this paper we study the transverse mass spectra and the HBT correlation 
functions of unidentified negatively charged hadrons (h~), which consist mainly of 
negative pions. 

2. The NA57 set-up and the data sample 

Descriptions of the NA57 apparatus can be found in references 0, [7J [8], and in 
reference [3] in particular for the 40 A GeV/c set-up. 

The tracking device of the NA57 experiment consisted of a telescope made of an 
array of silicon detector planes of 5x5 cm 2 cross-section placed in an approximately 
uniform magnetic field of 1.4 Tesla perpendicular to the beam line; the bulk of the 
detectors was closely packed in an approximately 30 cm long compact part used for 
pattern recognition. To improve the momentum resolution of high momentum tracks 
a lever arm detector (an array of four double-sided silicon micro-strip detectors) was 
placed downstream of the tracking telescope. 

The centrality of the Pb-Pb collisions is determined (off-line) by analyzing the 
charged particle multiplicity measured by two stations of micro-strip silicon detectors 
(MSD) which sample the pseudo-rapidity intervals 1.9 < 77 < 3 and 2.4 < 77 < 3.6. 

The results presented in this paper are based on the analysis of the data sample 
collected in Pb-Pb collisions at 40 A GeV/c. The selected sample of events corresponds 
to the most central 53% of the inelastic Pb-Pb cross-section and has been divided 
into five centrality classes (labelled with integers 0, 1, 2, 3, and 4, class 4 being the 
most central) according to the value of the charged particle multiplicity measured by 
the MSD. The procedure for the measurement of the multiplicity distribution and the 
determination of the collision centrality for each class is described in reference [9]. The 
fractions of the inelastic cross-section for the five classes, calculated assuming an inelastic 
Pb-Pb cross-section of 7.26 barn, are the same as those used for the study of strange 
particles [3] and are given in table [H 

Negatively charged tracks have been selected by requiring them to have clusters in 
more than 80% of the telescope planes and less than 30% of the clusters shared with 
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Table 1. Centrality ranges for the five classes. 



Class 


1 2 


3 


4 


a/o-ind (%) 40 to 53 


23 to 40 11 to 23 


4.5 to 11 


to 4.5 



other tracks, and using an impact parameter cut|§| to ensure they come from the main 
interaction vertex. 

3. Single particles m T spectra 

The acceptance region in the transverse momentum (p T ) versus rapidity (y) plane is 
shown in figure Dt^]. The limits of this window have been defined in order to exclude 
from the final sample the particles whose lines of flight are very close to the border 
limits of the telescope, where the systematic errors are more difficult to evaluate. 
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Figure 1. The y-pi acceptance window of negative hadrons. Dashed lines show the 
position of mid-rapidity (y cm = 2.225). 



A weight is assigned to each reconstructed h~ to correct for acceptance and 
reconstruction inefficiencies. The computational algorithm of the event weight is the 
same as that used for strange particles 0, [7|, [3]: a number of Monte Carlo events are 
generated, each event consisting of one simulated particle, with the same momentum 
of the real particle, merged with a real event of similar telescope hit multiplicity as the 
original event, and they are reconstructed with the same analysis tools as for real events. 
A total of about 10,000 h~, sampled uniformly over the full data taking periods, have 
been individually weighted with this method. 

§ The impact parameter is approximated as the distance from the primary vertex of the intersection 
of the measured particle trajectory with a plane transverse to the beam line passing through the target 
position. 

|| The rapidity has been evaluated assuming the pion mass. 



Expansion dynamics of Pb-Pb collisions at 40 A GeV/c 



5 



In order to check the stability of the results, the selection criteria (i.e. the impact 
parameter cut and the number of clusters associated to the track or shared with another 
track) have been varied, by changing their values one at a time. As a result of these 
studies we can estimate the contribution of the selection and correction procedure to 
the systematic errors on the slope of the l/m T diV/dm T distributions of h~ to be about 
7%. " 

The experimental procedure for the determination of the m T distribution is 
described in detail in references 0[3], where the results for strange particles are shown. 

The distribution of the double differential invariant cross-section — f 2N , has been 

assumed to factorize into a y and an m T dependent part: 

1 d 2 iV 1 diV 

— - A -7-{y,m T ) = f{y) — (m T ). (1) 

m T dm T ay m T dm T 

This assumption has been verified by considering the m T and y distributions for 

different slices, respectively, in rapidity and transverse mass. The shape of the rapidity 

distribution (f(y) in equation [I]) has been found to be well described by a Gaussian 

within our limited range. The hypotheses on the factorization of the double differential 

invariant cross-section (equation [1]) and on the shape of the rapidity distributions {f{y)) 

can introduce a contribution to the systematic error on the slope of the — 4^-(m T ) 

J 1 mi Ami v 1 ' 

distributions which has been estimated to be about 5%. 

The m T distribution has been parameterized as oc exp ^— yj- The i nverse 

slope parameter T app ("apparent temperature") has been extracted by means of a 
maximum likelihood fit of the measured double differential invariant cross-section to 
the formula 

1 d 2 N . . , / m T 



m T dm T dy V GXP \ T app J 
obtaining T app = 186 ± 2(stat) ± 16(syst) MeV for the most central 53% of the inelastic 
Pb-Pb cross-section. The differential invariant cross-section distribution is shown in 
figure [2] as a function of m T with superimposed the likelihood fit result. 

In the hydro-dynamical view, the apparent temperature is interpreted as due to the 
thermal motion coupled with a collective transverse flow of the fireball components [1], 
and it depends on mj^j]. At a given m T0 value, it can be calculated according to the 
formula pE]: 



2~app( m To) 



d dN . 
hm - — (bg— ^) 
mi-tniTo dm T arrij 



(3) 



This expression simplifies for m T ^> m: the apparent temperature is simply blue-shifted 
by the collective dynamics 



where T is the freeze-out temperature and (/3±) is the average transverse flow velocity. 

% In this view, the graphical interpretation of T app would be the inverse of the local tangent to the 
invariant l/m T diV/dm T distribution. See reference [7] for a detailed discussion. 
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Figure 2. Transverse mass spectrum of negatively charged particles for the most 
central 53% of the Pb-Pb inelastic cross-section. The superimposed exponential 
function has inverse slope equal to the T app value obtained from the maximum likelihood 
fit. 



3.1. Centrality dependence 

The transverse mass spectra of negative hadrons are shown in figure [3] for the five 
centrality classes defined in table [Q the inverse slope parameter T app does not depend 
on centrality and its values are given in table [2j 




Figure 3. Transverse mass spectra of negatively charged particles in Pb-Pb collisions 
at 40 A GeV/c for the five centrality classes of table [TJ The spectra of class 2, 3 and 
4 have been scaled by factors 2, 10 and 100, respectively, for display purposes. 
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Table 2. Inverse slopes (MeV) of the hit distributions of h~ in Pb-Pb collisions at 
40 A GeV/c for the five centrality classes defined in table [1] Systematic errors are 
estimated to be 8.5% for all centralities. 






1 


2 


3 


4 


186 ±6 


185 ±3 


185 ±3 


188 ±3 


187 ±4 



4. Two-particle correlation 

In this section we present the study of the Hanbury-Brown and Twiss correlation for 
negatively charged hadrons assumed to be pions (see section |4.1 .2p . An introduction to 
this topic and recent reviews of experimental results and theoretical developements can 
be found in references [TU| ITT] . 

The data will be analyzed in the framework of the hydro-dynamical inspired blast- 
wave model [TJ, which assumes cylindrical symmetry for an expanding fireball in local 
thermal equilibrium. Hydro-dynamics or parameterized models inspired by hydro- 
dynamics have shown to give a successful description of a number of observables, i.e. 
transverse momentum (p T ) and rapidity (y) distributions, direct and elliptical flow, two- 
particle correlation functions (for recent reviews see, e.g., references [12]). The single 
set of free parameters of the model to be obtained from the HBT study is: the kinetic 
freeze-out temperature (T), the average transverse flow velocity ((/3j_)), the Gaussian 
radius of the cylindrical system (Rg), the chaoticity of the pion emission (A), the proper 



time of the freeze-out (t/ = w tf — Xj) and the emission duration (At). 
Experimentally, the correlation function can be defined as 

where the signal S(q) is the measured distribution of the relative four-momentum 
q = pi — p 2 of two identical particles in one event and the background B(q) is the 
reference distribution built by pairing two particles taken from different events. For each 
pair of negative hadrons in the signal, there are about 15 pairs in the background, formed 
from events of similar multiplicity, so that the error on C2 is statistically dominated by 
the signal. The normalization factor iV is obtained by imposing that the integral of 
6*2(9) be equal to unit in the region of large q, where there is neither quantum (i.e. 
Bose-Einstein) nor other (e.g. Coulomb) correlations. 

The relative momentum q is measured on a pair-by-pair basis relative to the out- 
side-long reference, which is a Cartesian system defined by choosing the long axis along 
the beam direction; the out and side axes lay in the transverse plane with the former 
aligned with the average transverse momentum of the pair and the latter perpendicular 
to the other two axes. With respect to the laboratory system, the out-side-long reference 
is boosted pair-by-pair along the long axis in such a way to bring to rest the pair 
in that direction: in the literature this boosted system is usually referred to as the 
Longitudinally Co- Moving System (LCMS). In this reference system, the correlation 
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Figure 4. Two-particle acceptance window used in the Hanbury-Brown and Twiss 
analysis. The correlation function has been measured in each of the 20 rectangles 
patching the window together. The Kj and Y nn average values, measured over the 
rectangles, are reported in table [3J The dotted line is drawn at = y cm . 



function is parameterized according to a three-dimensional Gaussian function modified 
by the addition of an out-long cross-term: 

C 2 = 1 + A exp [-R 2 out q 2 - R 2 side q 2 s - i?L g ?i - 2\R 0l \R 0l q 0qi \ (6) 

Such a parameterization is suited for central collisions which are azimuthally symmetric. 
In case of a peripheral collision, the azimuthal symmetry is evidently broken; however 
it is recovered when building the correlation function from particles accumulated over 
many events with random impact parameters, which is our approach. The A parameter 
in equation [6] is referred to as the chaoticity parameter and it is expected to range in 
the interval [0, 1]. 



4-1. Study of the event sample properties 

4-1.1. Acceptance region The determination of the size and dynamical evolution of the 
system at freeze-out requires information about the dependence of the HBT radii on 
the mean momentum K = \{p\ + P2) of the pair. This can be parameterized by its 
transverse component and the pair rapidity (computed assuming the pion mass): 

K T = - J fjB.il + TV?) 2 + (Wl + Pz2) 2 (7) 



Y 



(Pyl +Py2f 

Ei + Ei 



— log 

2 & E 1 



f (Pzl 
Pxl + 



+ Pz2) 
Px2 



(8) 

E 2 ~ Pxl ~ Vx2 

With this in mind, our (two-particle) acceptance window, which is shown in figure HI 
has been binned in 20 rectangles in each of which the correlation function C^ifl) has 
been measured independently Table [3] displays the average values of the K T and Y nn 
distributions in each of the 20 bins used for the analysis. 

Additionally, in order to investigate the dynamics of the collision as a function of 
centrality, the whole analysis has been repeated in three centrality classes obtained from 
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Table 3. Averages of the Kj (uppermost number in the cells, in GeV/c units) and 
Y nn (bottom numbers, in blue on-line) distributions in the 20 rectangles of figure |4] 





1.157 


1.162 


1.163 


1.175 




2.65 


2.75 


2.85 


2.99 




0.804 


0.808 


0.809 


0.813 




2.64 


2.75 


2.85 


2.98 


T 


0.610 


0.615 


0.617 


0.620 


2.57 


2.68 


2.82 


2.97 


K T 


0.460 


0.468 


0.472 


0.478 


2.55 


2.68 


2.82 


2.97 




0.305 


0.332 


0.345 


0.358 




2.51 


2.67 


2.81 


2.95 




y > 

J 7T7T 7 



those defined in tabled) central (3 + 4), semi-central (2) and semi-peripheral (0 + 
which correspond to the ranges of centrality 0-11%, 11-23% and 23-53%, respectivel 

4-1-2. Contamination from non 7t~-iy~ pairs The correlation analysis presented here 
is based on pairs of negative hadrons (h~-h - ), which are expected to be dominated by 
pairs of identical pions. The purity decreases with increasing K T : based on the results 
of NA49 [13] we estimate that at K T w 0.3 GeV/c it is about 90%, at K T w 0.8 (1.2) 
GeV/c it goes down to about 80% (65%). The main contamination comes from ir~- 
K~pairs (from about 5% at K T pa 0.3 GeV/c to about 25% at K T pa 1.2 GeV/c) and also 
from K~-K~pairs (up to ~ 10% at K T pa 1.2 GeV/c). Other combinations, e.g. 7r~-p 
or K~-p, are negligible, the largest being about a few per cent for 7r~-p at K T pa 1.2 
GeV/c. 

Excepting K~-K~pairs at high K T , misidentified particles thus lead to counting 
of unlike pairs which do not give rise to Bose-Einstein correlations. Moreover, in Pb- 
Pb collisions at -Ebeam = 158 GeV per nucleon the measured HBT radius parameters 
of the K~-K~ (and K + -K + ) correlation functions [H] (NA49 Collaboration) are fully 
consistent with the published pion results and the hydro-dynamic expansion model. It 
has been shown in [15] (NA35 experiment), in [16] (NA49 experiment) and in [17] (WA97 
experiment) that the main effect of particle misidentification is to reduce the value of 
the chaoticity parameter A. Apart from particle misidentification which produces the 
strongest bias, the A parameter, shown in figure[5]as a function of K T , is also affected in a 
non-trivial way by resonance decays and other effects [IB] , which are expected to depend 
on the pair momentum. The A parameter, however, is not used for the reconstruction of 
the size and dynamical state of the source; the relevant source parameters are affected 
to the level of a few per cent only. 

+ For reasons of statistics the centrality classes and 1, as well as 3 and 4, have been combined. 




Expansion dynamics of Pb-Pb collisions at 40 A GeV/c 



10 



■ 0.8 
0.7 - 
0.6 - 
0.5 
0.4 

0.3 - 
0.2 - 



0.2 0.4 0.6 0.8 



K T ( GeV/c ) 



Figure 5. The chaoticity parameter A as a function of Kj for the most central 53% 
of the Pb-Pb cross-section. Statistical errors are shown. 



4-1.3. Two-track resolution The detector's ability to distinguish a pair of close tracks 
from a single track decreases with decreasing track separation, depending on the Si-pixel 
size. The correction for this inefficiency, applied pair-by-pair to the signal distribution 
S(q), was calculated via a Monte Carlo simulation by reconstructing a sample of 
(generated) two-tracks events embedded in real ones. The maximum inefficiency was 
found to be ~10% (relative to high q pairs) for pairs populating the 10-MeV-wide bin 
of the correlation functions centred at q = and it goes to zero outside that bin. 

4.1.4. Track splitting The possibility of ghost tracks in the data sample has been 
investigated. Fake tracks might be produced by the reconstruction algorithm due to, 
e.g., allowing track candidates to share several clusters. To avoid such a bias, tracks 
have been selected by requiring them to have clusters in more than 80% of the telescope 
planes and less than 30% of the clusters shared with other tracks^- As a final check, 
Monte Carlo events have been simulated by generating a sample of ir~ distributed in 
phase-space according to the measured double differential cross-section of equation [2j 
Then, the simulated events have been reconstructed using the same algorithm as for 
real events. Three tracks only have been reconstructed with one associated ghost out of 
250K generated tracks which pass through the pixel telescope. This effect is therefore 
negligible. 

4-1.5. Momentum resolution The quality of a correlation measurement is determined 
by the resolution of the two-particle momentum difference ("relative momentum") q. 
The relative momentum resolution was studied using a Monte Carlo chain based on 
the simulation code GEANT [19]. As listed in table HI we have estimated that the 
relative momentum projections q s and q 1 (the latter evaluated in the LCMS system) are 

* The same requirement was applied to the tracks used for the transverse mass spectra analysis 
discussed in section [3j as already mentioned. 
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Table 4. The resolution of relative momentum components q , q s and q± (MeV/c) in 
the LCMS system MS M function of K-y (GeV/c). 





K T < 0.5 


0.5 < K T < 1.0 


K T > 1.0 




22 


52 


110 


Qs 


15 


18 


21 


Qi 


10 


11 


13 



measured in this analysis with an error not larger than about 20MeV/c over our two- 
particle acceptance window. The q projection, on the other hand, shows a stronger K T 
dependence, which can introduce an important contribution to the systematics errors 
on the extracted HBT radii, as discussed below. 

4-1.6. Coulomb correction The observed two-particle correlation is expected to result 
from two different contributions, the Bose-Einstein effect and the Coulomb interaction. 
The Coulomb interaction between particles of same charge sign is repulsive, thus 
depleting the two-particle correlation function at small relative momenta. Two methods 
have been considered to correct for the Coulomb interaction. 

A standard procedure consists in applying g inv -dependent weights to each pair in 
the background distribution B(q) to get the Coulomb corrected correlation function: 

c ™ = "TOW = (9) 

1 + Aexp [-R 2 QUt q 2 - R 2 side ql - R\ ong ql - 2| j R o1 | j R o1 <?o<?i] 
where -ftTcoui is the squared Coulomb wave-function integrated over the whole source and 

Bowler and Sinyukov pointed out [20] that the correction procedure should be 
different in the case where not all particle pairs in the signal distribution are subject 
to the Coulomb correlation. The argument here is that if physics (or detector) related 
effects lead to a reduction of the observed quantum correlation strength, then the same 
effects also lead to a similar reduction of the Coulomb repulsion. Therefore, the strength 
of the Coulomb correction applied to the data should be linked to the experimentally 
observed A parameter. This method has been implemented recently by the CERES [2Tj, 
PHOBOS [22] and STAR [23] Collaborations. The correlation function in this procedure 
is fitted to 

C 2 (q) = ||| = (1 - A) + (10) 



A ■ i^coui(gmv) • [1 + exp(-R 2 ont ql - R 2 side q 2 - Rl ong q 2 - 2|i2 0l |i2 0l g g 1 ] 

where ifcoui(<2mv) is the same as in the standard procedure. We argue that this 
alternative method is appropriate when the correlation function is measured with a 
negligible contamination of non-identical particle pairs; in fact the particles of these 
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Figure 6. Coulomb correction for an expanding source (full triangles, see text for 
details) and for a point-like source (open squares) as compared to the Gamow factor 
(full line). 

pairs would Coulomb interact in the signal and would not in the background distribution. 
Therefore we consider it only as an alternative method to control possible systematic 
errors on the HBT radii due to the Coulomb correction, as discussed later. 

For both procedures, three kinds of pion emitting sources have been considered 
with a method similar to that described in reference [17] : 

• a point-like source (in space-time); 

• a static Gaussian source; 

• an expanding source, parameterized according to the same blast-wave model [I] 
used to describe the data. 

The point-like source was used for comparison with the analytical Gamow 
factor [2H [25] for the sake of validation of the method. A static Gaussian source with a 
radius of 5 fm yields a Coulomb correction smaller than the Gamow factor and slightly 
larger than that obtained for the expanding source. The computation of the Coulomb 
correction for the expanding source requires knowledge of the source parameters, which 
are not known a priori but should be obtained from the correlation study. Therefore 
an iterative procedure has been used, as discussed in reference [T7J, starting with the 
following values for the blast- wave model freeze-out parameters (to be discussed later): 
T = 120 MeV, (J3±) = 0.4, R G = 5 fm, r f = 5 fm/c and At = 1 fm/c. The results for 
the model parameters become stable after the first iteration^. 

In figure [6] the correction factor i^coui for an expanding source, which has been used 
in this analysis, is shown superimposed to the correction for a point-like source and to 
the Gamow factor. 

jj Varying the blast-wave model parameters by 20% produces a maximum variation of the Coulomb 
correction of about 5%. 
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4-2. Multidimensional fit 

Equation [6] has been fitted to the corrected correlation functions by maximizing the 
negative logarithmic likelihood function: 

-2 log L{R) = 2Yj[C l B l -S l \og{C l B i )+\og{S l \)} (11) 

i 

where Cj is the theoretical value of the correlation function for a given set of parameters 
R = (A, i? ou t, -Riong, -Rside? R i) m the i th bin of q and Si and Bi are, respectively, the 
distribution of signal and background in that bin. 

Two checks on the fit quality have been performed. First, the xV n df -values 
are calculated on the three-dimensional correlation function using the parameters 
corresponding to the maximum likelihood fit and are found to be distributed around 
unity. Second, the projections of the 3-dimensional correlation function onto each 
momentum difference component, with narrow cuts on the other components, have been 
fitted by the least-squares method to a Gaussian function, yielding results consistent 
with the results of the 3-dimensional fits. A typical sample of these projections is 
shown in figure [7] where the correlation function has been integrated along the other qj 
components in the intervals \qj\ < 30MeV. 
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Figure 7. Projections of the h~-h~ correlation function measured using Cartesian 
parameterization for pair rapidities 2.30 < < 3.25 and pair transverse momenta 
0.40 < Kj < 0.55 GeV/c. Full curves are the result of a least-squares fit to a Gaussian 
function (see the text). 



4-3. Systematic errors 

The main sources of systematic uncertainties on the HBT radii are attributed to three 
previously mentioned effects, namely, with decreasing importance, (i) the momentum 
resolution, (ii) the Coulomb correction procedure and (Hi) the double-track resolution. 

The systematic error due to the finite momentum resolution, which depends on 
the transverse momentum of the pair, has been evaluated by a Monte Carlo simulation 
with the following approach. A sample of pion pairs was generated, whose particles are 
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Table 5. Percentages of systematic errors on the HBT radii for three intervals of Kj. 





Kj < 0.5 


0.5 < Kj < 1.0 


Kj > 1.0 


Rout 


12 


14 


16 


Rside 


10 


10.5 


11 


R-long 


11 


11.5 


13 


Rol 


11.5 


15 


28 



extracted from the phase-space dplx dx „ of the expanding model with the parameters 
specified in section 14. 1 .61 The momentum- dependent part of this particle phase- 
space distribution is then propagated through the apparatus to take into account 
the experimental acceptance and reconstruction inefficency, in order to end up with 
the raw measured momentum distribution ■ Quantum and other correlations 

are introduced in a second stage using the method outlined in reference [26]. The 
same procedure is repeated for a second sample of pairs, obtained from the first one 
by smearing the momentum distribution according to the experimental momentum 
resolution. The variations of the extracted HBT radii provide an estimate of the 
systematic errors. 

As discussed in section 14.1.61 two Coulomb correction procedures have been 
implemented, each with different assumptions about the source properties (its space- 
time extent, static or expanding source, etc.). The differences in the HBT radii obtained 
with the different methods allow us to estimate a contribution to the systematic error 
of 10% for i? out , -Rside and -Riong, an d of 5% for the cross-term R 0l . 

The last source of systematical errors can be associated to the double-track 
resolution and arises from the finite statistics of Monte Carlo events generated to correct 
for this effect. It is of the order of a few per cent for all the radii. 

In table [5] the total systematic errors are given in percentages as a function of K T . 

5. Results and discussion 

5.1. Longitudinal expansion 

The cross-term R 0l provides information about the longitudinal expansion of the system; 
e.g. for a longitudinal Bjorken-like expansion [27] -R 0l should vanish in the LCMS 
system at mid-rapidity [28] . Figure M shows the K T dependence of the _R i parameters 
for the most central 53% of the Pb-Pb inelastic cross-section ("all") and for the three 
centrality classes specified at the end of section 14.1.11 The statistical and systematic 
errors are indicated in these and following plots with bars and shadow boxes (green 
on-line), respectively. Unfortunately, in most of the recent experimental works about 
HBT correlations in heavy ion collisions scarce consideration is given to the study of this 
parameter. Indeed, the cross term is often not used in the Cartesian parameterization 
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Figure 8. The R a i cross-term as a function of A' T for the most central 53% of the Pb- 
Pb cross-section (top- left panel) and for individual centrality classes 0-1 (top-right). 
2 (bottom-left) and 3-4 (bottom right) of tabic 1. The bars and the boxes show the 
statistical and systematic errors, respectively. 



(equation E]) or, when used, it is frequently assumed to be of a given sign (e.g. positive 
in the LCMS system). Such an approach could be justified at RHIC energies near mid- 
rapidity, where the Bjorken model should provide a good description of the longitudinal 
expansion; it is not at lower energy. Instead, the longitudinal dynamics is derived from 
the Yano-Koonin-Podgoretskii parameterization [29] of the HBT correlation function 



C 2 (q) 



Aexp [-R\ql_ - 7y k (9|| - ^yk<?o) 2 #j| - 7y k (<?o - v Yk q\\) 2 R 2 ] 



(12) 



where R±, R\\, R are called the "YKP radii'fl although R has temporal dimension 
(fm/c); -Uyk, called "Yano-Koonin velocity", is measured in units of c and 7 y k = 
'1 — td,)" 1 / 2 . In fact, in this parameterization the v yk parameter can be easily 



related, under some approximations, to the longitudinal flow velocity of the expanding 
source [301 EI]. 

The set of parameters of the Yano-Koonin parameterization can be computed from 
that of the Cartesian parameterization (and vice versa) via analytical relations [301 l31j 



ff The subscripts _L and || stand for "perpendicular" and "parallel" to the beam direction, respectively. 
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. In particular, for the Yano-Koonin velocity the relation is: 
A + B 



Vk 



2C 



2C 



A + B 



(13) 



where in the LCMS system, due to the vanishing of the longitudinal pair velocity, 
A = ° nt side , B = -Ri ong and C = — |-R i|-Roi; finally /?t is the transverse component of 

the pair velocity (3 ~ K/K°. 

From the Yano-Koonin velocity, one can compute its associated rapidity Y y k = 
I ln( 1 ^ y ^ ), which is often plotted as a function of the pair rapidity to study the 
longitudinal expansion (see, e.g., references [161 [22]). 111 figure [9] we plot this quantity, 
evaluated in the laboratory system, as a function of the pair rapidity Y nn also given in 
the same frame: Y y ^ b = Y^ MS + Y n7T . Statistical and systematic errors, shown at the 
la - confidence level, have been propagated from those of the out-side-long HBT radii. In 
these plots the correlation functions have been integrated in the range 0.1 < K T < 0.7 
GeV/c where the (K T ) has been evaluated to be 0.36, 0.41, 0.51 and 0.58 GeV/c, for 
the data points at Y n7T = 2.53, 2.68, 2.82 and 2.97, respectively. 

The source rapidity scales with the rapidity of the pair, indicating the presence of 
strong position-momentum correlations. A static source would exhibit no correlation 
and would correspond to a horizontal line at the rapidity of the center of mass. 
(Y^ k = y cm ). A source with strong dynamical correlations would correspond to a straight 
line along Y jk = Y n7T . The data are consistent with the latter scenario: particles emitted 
at a given rapidity are produced by a source moving collectively at the same rapidity. 



5.2. Transverse expansion 

The independence of i? S ide provides information about the transverse expansion of the 
source and its geometrical transverse size at freeze-out |33j . 

In figure [10] we plot -R s ide as a function of K T for the integrated centrality range 
and for the three centrality classes used in this HBT analysis. In the hydro-dynamical 
view, the decrease of i? S ide with increasing K T is due to the collective expansion in 
the transverse direction; the dependence of i? S ide on the pair momentum K can be 
parameterized by the formula [33l [M] 

J Rside(^T, = 1 (14) 

yi + MT^cosh^k-y^) 

where Mr = \/ Kj + m^, Rq is equal to the transverse geometric (Gaussian) radius of 
the source times l/v^2, T is the freeze-out temperature and (3± the slope of the (linear) 
transverse flow velocity profile: (3±(r) = fl±j^- A fit of equation [TU to the experimental 

data points provides the model parameter Rq and the ratio they are given in table [6] 
The measurement of the ratio -jr determines an allowed region for the pair of 
variables (/3±) and T, the former being computed from the linear slope /3± by assuming 
a uniform particle density of the source. Those allowed regions, at the la confidence 
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Figure 9. The Yano-Koonin rapidity Yyk as a function of the pair rapidity Y„ n 
measured in the interval 0.1 < Kj < 0.7 GeV/c (see text for the averages of Kj in 
each bin) . The most central 53% of the Pb-Pb cross-section (top- left panel) and the 
individual ccntrality classes 0-1 (top-right), 2 (bottom-left) and 3-4 (bottom right) of 
table 1 are shown. Full circles are data, open circles are data reflected about mid- 
rapidity. Statistical and systematic errors are indicated with line bars and shaded 
boxes (green on-line), respectively. 

Table 6. The transverse geometric parameter Rq and the ratio 0]_/T, T being the 
freeze-out temperature. 



Centrality 


Rg (fin) 


^ (GeV- 1 ) 


ALL 


5.7 ±0.9 


2.9 ± 1.4 


0-1 


4.3 ±0.9 


2.7 ± 1.7 


2 


6.9 ±2.2 


5.4 ±3.3 


3-4 


6.6 ± 1.5 


4.1 ±2.6 



level, are shown in figure [IT] they correspond to the wide bands with positive slopes (in 
black on-line). 
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Figure 10. The R S id e radius as a function of K T for the most central 53% of the Pb- 
Pb cross-section (top-left panel) and for individual centrality classes 0-1 (top- right), 2 
(bottom- left) and 3-4 (bottom right) of table 1. Statistical and systematic errors arc 
indicated with line bars and shaded boxes (green on-line), respectively. 



In order to disentangle the two parameters, we followed a technique first used by the 
NA49 Collaboration [16], where the measured single-particle m T spectra are exploited. 
In this approach, the inverse slope parameters of the m T distributions are interpreted 
according to the blue-shift formula HI as discussed in section [31 This independent 
measurement provides another allowed region in the freeze-out parameter space, which 
corresponds to the narrow bands with negative slopes of figure [Til (blue on-line). The 
disentangled values for T and ((3±) are reported in table [3 

In figure [TTJ we also show the la confidence regions obtained from the blast- wave 
analysis of the m T spectra of singly-strange particles (Kg, A and A) [3]. They correspond 
to the small closed contours (red on-line) of figure [UJ with the markers indicating the 
optimal fit locations. Our results thus suggest compatible freeze-out conditions for 
singly-strange particles and h~ (mainly negative pions) in Pb-Pb collisions at 40 A 
GeV/c. On the other hand, the analysis of the transverse mass spectra of multiply- 
strange hyperons (S - and Q~) suggested that for the same colliding system these 
particles may undergo an earlier freeze-out than singly-strange particles [3]. 
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Figure 11. Contour plots at la confidence level (statistical error) in the (f3±)-T 
freeze-out parameter space, see text for details. 



Table 7. The freeze-out temperature T and the average transverse flow velocity (/3j_) 
as obtained from the study of negatively charged hadrons. 



Centrality 


T(MeV) 


(P±) 


ALL 
0-1 

2 

3-4 


122^ 
114^ 


0A0t° % 
0.371°;}°, 
0.501°;}° 
0.45l°i° 



Postponing to the next section a global discussion about these results, we note here 
that the most peripheral class 0-1 features a smaller transverse freeze-out radius, a lower 
transverse expansion velocity and a higher temperature. This suggests an expansion 
on a smaller scale; the higher temperature at the freeze-out may be interpreted as the 
remnant of an earlier decoupling of the expanding system. Such a centrality dependence 
of the freeze-out parameters is well established at the SPS and RHIC both from studies 
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of the transverse mass spectra [3J El ESJ ESI EZ] and from HBT analyses [T7J [2TJ |2"3"1 [55] . 

5. 3. Temporal characterization of the expansion 

Two parameters characterize the temporal evolution of the expansion dynamics, the 
proper time r of the kinetic freeze-out and the duration of the pion emission At. 



5.3.1. Proper time of freeze-out Information about the evolution time-scale of the 
source, or proper time of freeze-out, can be extracted from the independence of the 
i? long radius. In figure [12] we show the K T dependence of R long . We have fitted this 
parameter to a formula first suggested by Sinyukov and collaborators [39] and then 
improved by Chapman et al [33] : 

iW^x, Y™) = \P~ h( J { v \ (15) 
V m T cosh(F yk - Ym) 

This formula assumes an instantaneous freeze-out in proper time (i.e., Ar = 0); such an 
approximation is justified by the small Ar found in the present (see later) and similar 
analyses. 




Figure 12. The Ri ong radius as a function of K T . Statistical and systematic errors 
are indicated with line bars and shaded boxes (green on-line), respectively. 
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Table 8. The proper time of freeze-out (r) in integrated centrality class (all) and for 
individual centrality classes. The second error is associated to the uncertainties on the 
value of T, as discussed in the text. 



Centrality 


r (fm/c) 


ALL 


6.7 ±0.2 ±g;| 


0-1 


5.0 ±0.3 tol 


2 


7.2 ± 0.2 t° 5 7 


3-4 


7.5 ± 0.3 t° i 



The results of the bidimensional fits of equation [15] to the experimental data are 
presented in table [SJ In the fitting procedure we have used the freeze-out temperature 
determined from the analysis of the transverse dynamics, as discussed in section 15.21 
This introduces an error which is kept separate in table [H] from that due to the -Ri ong 
uncertainties. As observed for the transverse expansion, the class 0-1 shows an expansion 
on a smaller scale, its freeze-out (proper) time r being significantly shorter than for the 
two most central classes. 

5.3.2. Mean duration of the pion emission It has been proposed that the existence of a 
strong first order phase transition and an accordingly long-lived mixed phase would be 
observable by a large outward radius R out compared to -R si d e) indicating a long duration 
of the pion emission At 001 El 021 US M, HH] : 

At 2 = ° ut R2 side (16) 
Pt 

where /3t is the average transverse velocity of the pair. The dependence of the R out 
parameter on K T is shown in figure [T3"l We found that the ratio -R ut/-Rside is compatible 
with one in the explored K T range (figure [T4"j) . Values smaller than unity for this ratio 
can be expected for sources with surface dominated emission |36] , such as emission from 
an expanding shell. We plot the quantity s^s f or the individual centrality classes 

in figure [T5l The data do not support the scenario of a long-lived source, but more likely 
that of a sudden freeze-out, At being of the order of a few fm/c for small K T and nearly 
zero for high K T . 

5.4. Consistency 

In the first three rows of table [9] we have summarized the model parameters obtained 
from the study of h~ single-particle m T spectra and two-particle correlation functions 
for the three centrality classes considered. These data provide implicitly a dynamical 
picture of the collision process. 

To see if this picture is self-consistent, we can compare the two-dimensional rms 
width i?^g ZCOUt = \/2Rg with the two-dimensional rms width of a cold lead nucleus 
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Figure 13. The Rout radius as a function of Kj for the most central 53% of the Pb- 
Pb cross-section (top-left panel) and for individual centrality classes 0-1 (top-right), 2 
(bottom- left) and 3-4 (bottom right) of table 1. Statistical and systematic errors are 
indicated with line bars and shaded boxes (green on-line), respectively. 

i?jjj s = y^l^A 1 / 3 ~ 4.5 fm. For the most central collisions (class 3-4), the system 
expands by a factor «2 or, equivalently, by about 5 fm in the transverse direction. If 
the transverse flow velocity is equal to (3± = §(/3j_) — 0.67 at the surface during the 
whole expansion, in a time of r ~ 7.5 fm/c nuclear matter would travel over Tq/3± ~ 5 
fm in the transverse direction. This is consistent with the previous estimate from the 
difference i?^g Zeout — i?^ s . Therefore, the dynamical description of the system expansion 
is internally consistent. Similar pictures can be drawn for the other two classes. 

6. Comparison with other experiments 

In figure [16] we show HBT radii versus K T measured in central Pb-Pb collisions at 40 
A GeV/c by various experiments. The overall agreement of the data is good for the 
Rside an d -Riong radii. A discrepancy is observed for the R out radius, with the CERES 
results significantly smaller than the others. The values of the parameter R o1 depend 
strongly on the pair rapidity, hence a comparison with the CERES data, which are given 
at slightly backward rapidity, would not be straightforward. On the other hand, the 
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Figure 14. The ratio -R ut /-Rside as a function of Kj for the data sample corresponding 
to the most central 53% of the Pb-Pb inelastic cross-section. 
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Figure 15. The At 2 parameter of equation [16] as a function of Kj for the three 
centrality classes. Statistical and systematic errors are indicated with line bars and 
shaded boxes (green on-line), respectively. 



agreement between NA49 and NA57 data, both shown at slightly forward rapidity over 
a similar rapidity range, is satisfactory also for the R o1 parameter. NA57 is the sole 
experiment that explores the high K r region and our data suggest that the HBT radii 
-Rside, -Rout and -Riong decrease steadily up to K T « 1.2 GeV/c. 

Superimposed on the experimental data, we show in figure [16] (with asterisks, 
connected by lines) results by Qingfeng Li et al. [H] based on the UrQMD v2.2 transport 
model [30] • The radii -Ri on g and -R s ide are reasonably well in line with experimental 
data; the predicted R out values are larger than those of the experimental data. As 
a consequence, the extracted quantity a/ B% ut — Rl ide of the pion emission source (see 
equation [TBI) becomes larger than the experimental evaluation. 

In figure [17] we show the plot of the Yano-Koonin rapidity versus the pair rapidity, 
both evaluated in the centre-of-mass rest frame, as measured at different energies. 
Figure [17] reveals a roughly universal dependence of Y Y ^ on Y nn for pions from central 
collisions, depending weakly, if at all, on ^/snn- This trend is particularly striking given 
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Coll. 


V^NN 


Centrality 


Rq (fm) 


T (MeV) 


(AO 


t (fm/c) 


At (fm/c) 


NA57 


8.8 


- 


11 % 


6.6 ± 1.5 


114li4 


0.451^ 


7.5 ± 0.7 




NA57 


8.8 


11 


- 23 % 


6.9 ± 2.2 


IO6I13 


0.50tg; 2 


7.2 ± 0.6 


1-2 


NA57 


8.8 


23 


- 53 % 


4.3 ± 0.9 


- ^ - -1-90 

1251?° 


0.37l°; 2 


5.0 ± 0.5 




CERES 


8.8 


- 


5 % 


„ rt w J_7 7 
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8.8 


5 - 


10 % 


y.y_ 2 


120 
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„ „_L1 1 
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_ rA _L(] Of) 
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in P. 

c n+u.t) 

o.y_o.4 
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0.50l°;°? 
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12.3 


15 
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17.3 


- 
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1 A 

6 9+ - 3 
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1 A A£ 
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8.2 ±0.1 
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5 - 


10 % 


3 

r n + U.O 

o.y_ 2 
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^ 1- r.4-0 On 


7.9 ± 0.1 


2-3 


CERES 


17.3 


10 


- 15 % 


5-9Io.4 
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17.3 


15 


- 30 % 


_ r+0.4 

O.0_ 4 


120 


„ £.£.+0.03 
U - 00 -0.03 


7.3 ± 0.1 




WA97 


17.3 


- 


5% 


5.1 ± 0.6 


1201 


„ . „4_() 07 
4fv 
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5.6 ± 0.2 




WA97 


17.3 


5 - 


12% 
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117l{? 


^ . ok 

0.48Io;n 
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w 
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17.3 


12 


- 25% 


4.6 ± 0.4 


-1 -4-1 ^ 

1211^ 


0.47tHo 
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WA97 


17.3 


25 


-40% 


3.2 ± 0.3 


140li3 


— ^ ~ 4-fl HQ 

u - ou -0.16 


3.7 ± 0.2 




NA49 


17.3 


- 


3% 


6.5 ± 0.5 


120 ± 12 


0.55 ±0.12 


8 


3-4 


T~\ TTT1A T~T~"\7~ 

PHENIX 


130 


- 


30% 


8.1 ± 0.3 


125 


0.51 






STAR 


200 


- 


5% 


13.3 ±0.2 


97 ±2 


0.59 ±0.05 


9.0 ±0.3 


2.83 ±0.19 


STAR 


200 


5 - 


10% 


12.6 ±0.2 


98 ±2 


0.58 ±0.05 


8.7 ±0.2 


2.45 ±0.17 


STAR 


200 


10 


- 20% 


11.5 ±0.2 


98 ±3 


0.56 ±0.06 


8.1 ±0.2 


2.35 ±0.16 


STAR 


200 


20 


- 30% 


10.5 ±0.1 


100 ±2 


0.54 ±0.06 


7.2 ±0.1 


2.10 ±0.09 


STAR 


200 


30 


- 50% 


8.8 ±0.1 


108 ±2 


0.51 ±0.06 


5.9 ±0.1 


1.74 ±0.12 


STAR 


200 


50 


- 80% 


6.5 ±0.1 


113 ±2 


0.41 ±0.06 


4.0 ±0.2 


1.73 ±0.10 



Table 9. Freeze-out parameters from HBT measurements in heavy ion collisions 
at SPS and RHIC. In the CERES [21] and PHENIX [38] analyses the freeze-out 
temperature has been fixed to the quoted values. The reported values of (/3j_) for the 
PHENIX results have been computed from their parameterization of the transverse 
flow, which is assumed to be linear in the rapidity profile. In the STAR analysis [23], 
the values of (/3±) and T have been extracted from blast-wave fits to pion, kaon, and 
proton transverse momentum spectra [54j and used to describe HBT radii. 

the very different centre-of-mass projectile rapidities (« 1.57 and 5.5 for a/snn = 5 and 
200 GeV, respectively) and the corresponding widths of the pion distributions dN/dy. 
The centrality dependence of HBT radii is shown in figure [18] for a wide range of 

collision energies. The left panels show the dependence on the number of participating 

1/3 

nucleons, N pait . All of the radii exhibit a linear scaling in N pait . Only the slope of the 
-Riong dependence shows a significant increase from the AGS to RHIC, consistent with 
a lifetime that increases with both centrality and ^/snn- The trend of increasing i?i on g 
with increasing ^/snn is reversed for ^/snn < 5 GeV [52] . We note that R out radii from 
CERES at 40 A GeV/c are well below the observed systematics. 

The right panels of figure [TBI show the same radii as a function of (dN^/dr)) 1 / 3 . 



Expansion dynamics of Pb-Pb collisions at 40 A GeV/c 



25 



{ 



J I CERES 
5 NA49 
j NAS7 
I UrQMD 



0.2 0.4 0.6 0.8 1 1.2 

Kt (GeV/c) 



j $ CERES 
} NA49 
. NAS7 
. UrQMD 



0.2 0.4 0.6 0.8 1 1.2 

K T (GeV/o) 



: □ 


J } CERES 




i NA49 






; NAS7 




- \ 


I UrQMD 








: i 




















I , , , 1 , , , ! , , , 1 , , 





0.2 0.4 0.6 0.8 1 1.2 

K T (GeV/c) 



- 
-1 - 



5 NA49 
. NA57 
; UrQMb 



4 



0.2 0.4 0.6 0.8 



K T (GeV/c) 



Figure 16. The HBT radii as a function of Kj in Pb-Pb collisions at 40 A GeV/c 
beam momentum. The CERES data set [H] corresponds to the most central 5% of the 
inelastic Au-Pb cross-section; those of NA49 [47] and NA57 correspond to the most 
central 7.2% and 11% of the inelastic Pb-Pb cross-section, respectively. 



The primary motivation for exploring the (dA^h/dr/) 1 / 3 dependence is its relation to 
the final state geometry through the density at freeze-out. However, the two scaling 
quantities are highly correlated. In fact, the values of dN^/dr] shown on the right 
side of figure [18] are derived (as suggested in reference [10]) from Ap art using the A"p art 
parameterizations given in [53], and conversely, the Ap art values are often calculated 
from multiplicity distributions using a Glauber model. Given this caveat, the i? S ide and 
i?ion g values exhibit a linear dependence on (diV^/d?]) 1 / 3 . The similar behaviour from 
v /snn of 5 to 200 GeV leads one to believe that the approximate N pait scaling (initial 
overlap geometry) is a result of the scaling with multiplicity (final freeze-out geometry) 
and not the other way around. 

In table [9] we present, along with our results, a compilation of results on the source 
parameters obtained at the SPS and RHIC from blast-wave fits to HBT correlation 
functions. At all collision energies it is observed that, with increasing centrality: (i) the 
geometrical transverse dimension of the system at freeze-out, the average transverse flow, 
the proper time of the freeze-out and the duration of pion emission increase; (ii) the 
thermal freeze-out temperature decreases. 

With increasing ^/snn, both the proper time of the freeze-out (which can be 
assumed as the duration of the expansion) and the average transverse flow velocity 
increase, thus resulting in a stronger increase of the transverse radius parameter Rq at 
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Figure 17. The Yano-Koonin rapidity plotted as a function of the pair rapidity for 
central Pb(Au)-Pb(Au) over a broad range of energies (open symbols are for tt~-tt~, 
closed symbols for ir + -ir + ). Both quantities are in the center of mass frame of the 
colliding system. Data of other experiments are taken from references [50l [16j H3 E2] ■ 

the freeze-out: after a faster expansion which lasts longer, the system ends up showing 
much larger spatial extent. The freeze-out temperature at top RHIC energy is also 
significantly smaller than that measured at SPS for a given collision centrality. This 
can be explained again as due to a longer and faster expansion, which causes the system 
to cool down more effectively. On the other hand, similar freeze-out emission durations 
of a few fm/c are observed independently of the collision energy. 

7. Conclusions 

The analyses of transverse mass spectra of negatively charged hadrons and of Hanbury- 
Brown and Twiss correlation functions have provided a description of the later stages of 
the expansion dynamics of the system formed in Pb-Pb collisions at 40 A GeV/c. They 
have been performed over a centrality range corresponding to the most central 53% of 
the inelastic Pb-Pb cross section. 

Based on a model derived from hydro-dynamics, we have determined parameters 
describing the spatial dimension and the dynamical/thermal state of the system, namely 
the transverse radius and the temperature of the fireball at kinetic freeze-out, the average 
velocity of the collective flow in the transverse direction, the proper time of the freeze-out 
and the mean duration of particle emission. For central collisions, the system undergoes 
a strong collective expansion in the transverse directions, with an average velocity of 
about 45% of the speed of light, which lasts for about 7 fm/c. After this expansion, the 
system has doubled its transverse size and it freezes out at a temperature of about 110 
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Figure 18. The HBT radius dependence on number of participants (left) and on 
charged particle multiplicity (right). Data are for Pb(Au)-Pb(Au) collisions at several 
values of ^/snn, and also for Si- A collisions at the lowest energy. Average transverse 
momentum (K T ) is ~450 MeV/c for the PHENIX data [38] and ~390 MeV/c for the 
others j5Tl [2T1 1X71 [23] . Open symbols are for 7r~-7r~, closed symbols for 7r + -7r + , but 
for the STAR, PHENIX and CERES results which combine the two sets. Systematic 
and statistical errors have been added when both available. 



MeV. For less central collisions, the expansion proceeds on a smaller scale; the higher 
temperature at freeze-out may be interpreted as the remnant of an earlier decoupling 
of the expanding system. 

Concerning the time evolution, we have deduced that h~ emission is a fast process 
(At ~ 1 — 2 fm/c) starting after 7.5 fm/c (5 fm/c) in the case of central (semi-peripheral) 
collisions, that takes places in the bulk (i.e. it is not surface dominated). This resembles 
the decoupling process of photons in the early universe. 

In the longitudinal direction, the indications are that particles emitted at a given 
rapidity are produced by source elements moving collectively at the same rapidity. 

The kinetic freeze-out conditions for h~, namely temperature and transverse flow, 
are compatible with those obtained, within the same model, by studying the m T spectra 
of singly-strange particles (Kg, A and A) [3]. 

Quite universal kinetic freeze-out conditions are observed with increasing ^/snn 
from SPS to RHIC, with moderate increases (by less than 20%) of the transverse flow 
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velocity and of the duration of the expansion. The latter effect is consistent with the 
lower (by ~20%) freeze-out temperatures measured at top RHIC energy (^/snn = 200 
GeV); on the other hand, the combined effect of a faster and longer expansion causes 
the final transverse size at freeze-out to be significantly larger (about a factor of two for 
central collisions) at top RHIC energy than at SPS. 
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